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ABSTRACT 


In this paper, we study biharmonic Riemannian submersions 7 : M? x R > 
(N?,h) from a product manifold onto a surface and obtain some local charac- 
terizations of such biharmonic maps. Our results show that when the target 
surface is flat, then a proper biharmonic Riemannian submersion 7 : M?xR > 
(N?, h) is locally a projection of a special twisted product, and when the tar- 
get surface is non-flat, m is locally a special map between two warped product 
spaces with a warping function that solves a single ODE. As a by-product, we 
also prove that there is a unique proper biharmonic Riemannian submersion 


H? x R > R? given by the projection of a warped product. 


1. INTRODUCTION AND PRELIMINARIES 


A harmonic map is a map y : (M,g) > (N,h) between Riemannian man- 
ifolds whose tension filed vanishes identically, i.e., T(g) = Trace,Vdy = 0. A 
biharmonic map is one whose bitension filed solves the PDEs 


(1) Ta(p) := Trace,(V?°V? — Vem)T(Y) — Trace, R (dy, T(y))dy = 0, 


where R is the curvature operator of (N,h) defined by 
RN(X,Y)Z = [V¥, V¥]Z — VikyZ. 


Clearly, any harmonic map is a biharmonic map. A biharmonic map which is not 
harmonic is called a proper biharmonic map. 


The geometric study of biharmonic maps focuses on biharmonicity of maps 
with geometric interest like isometric immersions or Riemannian submersions, 
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and the geometry and topology of the spaces related to the existence of such 
geometrically special biharmonic maps. For example, in the study of biharmonic 
submanifolds (biharmonic isometric immersions), a fundamental problem is to 
classify biharmonic submanifolds in certain model spaces. Although many pro- 
gresses have been made the following conjectures remain open in general cases. 


Chen’s Conjecture ([8]): Any biharmonic submanifolds in Euclidean space 
is minimal. 


Conjecture (Balmus-Montaldo-Oniciuc [3]): A biharmonic submanifold in 
sphere has constant mean curvature; and any proper biharmonic hypersurface in 
py is an open part of S™ (Zz) or the generalized Clifford torus S?(Z) x Sz) 
with p+q = m, p # q. 


For more detailed information and some recent progress on biharmonic sub- 
manifolds see the recent book [19] and the vast references therein. 

Riemannian submersions are a dual concept of isometric immersions (i.e., sub- 
manifolds). The study of biharmonicity of Riemannian submersions was initiated 
in [16]. A useful tool of using the so-called integrability data to study biharmonic 
Riemannian submersion from a generic 3-manifold was introduced in [21]. This 
was later generalized to higher dimensions with one dimensional fibers in [1]. 
Complete classifications of Riemannian submersions from a 3-dimensional space 
forms and more general BCV spaces into a surface were obtained in [21, 22]. 


In this paper, we study biharmonic submersions from the product space M? xR, 


where M? is a general 2-dimensional manifold. A reason for the choice of the 


product space M? x R is that it includes many well-known model spaces, such as 
S?xR, H?xR, R, and twisted spaces R° x2» R = (RÌ, dy? + dz? + e??@") dx). 
We give some local characterizations of biharmonic (including harmonic) Rieman- 


nians submersions from M? x R onto a surface. These include: a Riemannian 


submersion m : M? x R > (N?,h) is harmonic if and only if it is locally the 
projection onto the first factor followed by a Riemannian covering map (Theorem 


2.2); A Riemannian submersion m : M? x R —> R° is proper biharmonic, then 


it is locally a projection of a special twisted product (Theorem 2.9), and if a 


Riemannian submersion 7 : M? x R > (N?,h) into a non-flat surface is proper 
biharmonic, then 7 is locally a special map between two warped product spaces 


with a warping function that solves a single ODE (Corollary 2.10). 
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2. HARMONIC AND BIHARMONIC RIEMANNIAN SUBMERSIONS FROM M? x R 


We will use the following useful local orthonormal frame in the paper. 


Lemma 2.1. For any point on M? x R, there is a neighborhood U x R and 
local coordinates (t,s,z) such that the product metric on M? x R takes the form 
eats) dt? + ds? +dz?. Furthermore, the orthonormal frame {E = e (639, By = 
s, E3 = 0,} satisfies 


(2) |E, E2] = fE, all other |E, E3] = 0, i,j — 1, 2,3, 
Vg E = —f Eo, Va, Eo = fE, all other Ve Ej = 0, 1, J = 1, 2,3, 


where f = qs, and the only possible nonzero components of the (Ricci) curvatures 
of M? x R are given by 


(3) Ric (E, E) = Ric( E>, E») = Riz12 = g( R(E, E>) Ez, E1) = — E(f) = fa 
which is the Gauss curvature of M? on U denoted by K? = —E,(f) — f?. 


Proof. It is well known that around any point of a 2-dimensional Riemannian 
manifold (M?, g) there exists a semi-geodesic coordinate system so that locally it 
can be described as (U C M?, g = e744) dt? + ds?). One can easily checked that 
{E, = e~%49)0,, Ey = 0,, Ez = 0,} is an orthonormal frame satisfies (2) and (3) 
with f = qs. 


Remark 1. Note that the Lemma 2.1 actually implies that any local semi-geodesic 
coordinate system gives rise to an orthonormal frame {F;} with the properties 
(2) and (3). For latter reference we call such an orthonormal frame semi-geodesic 
orthonormal frame. 


Let 7: M? xR > (N?,h) be a Riemannian submersion, and let {e1, e2, e3} be 
a local orthonormal frame on W C M?xR with eg vertical. Let {fi, fo, f3, K1, K2, O} 
be the generalized integrability associated to this frame. Note that f3 = 0 if and 


only if the above frame is adapted to 7. 
3 
Using the relation e; = J` alE;, as it was computed in [22], we have the 
j=l 


following Lie brackets and the components of the Levi-Civita connection with 
respect to this frame as 


(4) 


[e1, e3] = fse2 + K1e3, [e2,e3] = —fze1 + K2e3, [e1,€2] = frei + foes — 2063, 
Vee = — fiez, Ve e2 = fier — 0€3, Ve e3 = 0€2, 
Ver€1 = —foe2+oe3, Vero = frei, Verses = —9€1, 


Vex€1 = —Ki1€3 + (0 — f3)e2, Veze2 = — (0 — f3)e1 — K2€3, Vege3 = K1€1 + Kea, 
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and the only possible nonzero components of the Riemannian curvature R of 
M? x Ras 


R(e1, €3, €1, €2) = —€1(0) + 2kyo = —a3a3K™’, 
R(e1,€3, €1, €3) = e1(K1) + 0? — K? + kafi = ($P KË, 
R(e1, €3, €2, €3) = €1(K2) — e3(0) — ki fi — Kike = -ofa K”, 
(5) 4 R(e1,€2,€1,€2) = e1(f2) — e2(f1) — F? — J3 + 2fgo — 30? = (ag) KP, 
R(e1, €2, €2, €3) = —e2(0) + 2K20 = BK, 
( = 


€2(K1) + €3(0) + Kof — Kik = -03a K, 


R(e2, €3, €2, €3) = o? + €2(K2) — Kı f2 — K3 = (RE, 


Ria, €3, €1, €3 


Finally, the Gauss curvature of the base space is given by 
(6) K™ = e1(fo) — ea(fi) — fi — £3 + 2fao. 


Remark 2. From (4) we see that øo = 0 on the neighborhood W if and only if the 
horizontal distribution of the Riemannian submersion 7 on W is integrable. 


2.1. Harmonic Riemannian submersions from M?xR. Harmonic Riemann- 


ian submersions M? — N? are a special subclass of horizontally homothetic har- 
monic morphisms with totally geodesic fibers. For general harmonic morphisms 
and their applications and interesting links to other areas of mathematics see the 
book [4]. It is well known that the only Riemannian submersion from $° is the 
Hopf fibration S3 — S$? which is harmonic since the fibers are totally geodes- 
ic. Using the Bernstein theorem for harmonic morphism (see Baird-Wood [5]), 
one can easily deduce (see also [21]) that there is no harmonic Riemannian sub- 
mersion H? — (N?,h) from a 3-dimensional hyperbolic space form, and that any 
globally defined harmonic Riemannian submersion ¢ : R? — (N?, h) is an orthog- 
onal projection R? > R° followed by a Riemannian covering map R? > (N?, h). 
For some recent work on the classifications of harmonic Riemannian submersions 
from Thurston’s 3-dimensional geometries, BCV 3-spaces, and Berger 3-spheres 
see [23]. In this subsection, we study harmonic Riemannian submersions from 
M? xR. 


Theorem 2.2. A Riemannian submersion 7: M? x R > (N?,h) is harmonic if 
and only if locally it is, up to an isometry of the domain, the projection onto the 


first factor followed by a Riemannian covering map. 


Proof. By Lemma 2.1, for any point on M? x R, there is a neighborhood W = 


3o 
U x R on which we have an orthonormal frame {E;}. Let e; = X, al E; be a 
j=1 


local orthonormal frame with e3 vertical. Using Proposition 2.2 in [23] and (5), 
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we conclude that if the Riemannian submersion 7 : M? x RD W = (N?,h) is 
harmonic, then we have 


(7) 0? =(a8P?K™ = (a2)?K™’, abab KM = 0, and KN = 30? + (a3)?K™. 


By continuity of KM’, we may assume that either KY? = 0 or KM? Æ 0 on 
a neighborhood denoted by W by an abuse of notation. For the first case, (7) 
implies that o = 0, and KY = KM? = 0 on W. For the second case, we use 
K™ Æ 0 on W together with the 1st and the 2nd equation of (7) to conclude 
that a? = a3 = 0 and hence (a3)? = 1, which also implies that o = 0 and hence 
the Gauss curvature of the base space KN = o? + (a3)?K™ = (a3)?K™ = KW. 
Thus, in either case, we have o = 0 on a neighborhood W which means the 
Riemannian submersion m : M? x R D W —> (N?,h) has integrable horizon- 
tal distribution. This, together with the assumption that m|w is harmonic and 


hence it has totally geodesic fibers, implies that m|w is a Riemannian submersion 
with totally geodesic fibers and integrable horizontal distribution whose inte- 
gral submanifold is an open neighborhood U of M? that is isometric to an open 
neighborhood of (N?, h). Therefore, t|w is the projection along the fibers onto 


U followed by a Riemannian covering map onto V C N?. 


2.2. Biharmonic Riemannian submersions from M? xR. In this subsection, 


we will characterize all proper biharmonic Riemannian submersions from M? x R, 
which are not harmonic. 
We will use the following lemma in the rest of the paper. 


Lemma 2.3. ([21]) Let 7: (M?,g) —> (N?,h) be a Riemannian submersion with 
an adapted frame {e1, e2, e3} and the integrability data fi, fo,K1, K2 and ø. 
Then, the Riemannian submersion m is biharmonic if and only if 


(8) 
-AMK -2D fuel) = ro Ds lelh) — kifi) + K1 (-x* + 2 =0 


i=l 


—A™ ko +2 3 fiei(K1) + Ky Tlel) = Ki fi) + Ko (-« + =f] — 0, 


where KN = RN. Om = €1(fo) — eo( f1) — f? — f2 is the Gauss curvature of 
(N2, h). 


Lemma 2.4. (see [22]) Let 7 : (M?, g) > (N?, h) be a Riemannian submersion 
from Riemannian 3-manifolds and {e1, e2,e3} be any local orthonormal frame with 
e3 tangent to the fibers. If Ve,e1 = 0, then either V-,e2 = 0; or V-,e2 £ 0, and 
the frame {e1, €2,e3} is adapted to the Riemannian submersion r. 
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The main tool used to prove our main theorems is a special orthonormal frame 
adapted to the Riemannian submersion. First we prove the existence of such 
frame. 


Proposition 2.5. Let 7: M? x R > (N?,h) be a Riemannian submersion and 
{Fy, Eo, E3 = 2} be a local semi-geodesic orthonormal frame stated in Lemma 
2.1. Then, there exists an orthonormal frame 

e1 = cos 0 E1 + sin 0 E», 
(9) e = — cos Q (sin 0E; — cos 0 E3) + sin a E3, 

e3 = sina(sin 0E; — cos 0 E2) + cos a Es, 


such that e3 is vertical, Vee1 = 0, where cosd = (e1, E1), cosa = (e3, E3). 
Furthermore, the generalized integrability data of {e;} are given by 


fi=0, fo= — f cos? a — sina cos aF3(), fs = —E3(6), 


(10) o = —f sina cosa — sin? vE3(@), kı = — f sin? a + sina cos a E3(0), 
k2 = sin ae3(cosa@) — cos ae3(sina) = —e3(a), 

where 

(11) f = -sin 0E (0) + cos 0E2(0) + f sin 0. 


Proof. Let e3 be the unit vector field tangent to the fibers of m. Clearly, if 
e3 = +E3, then e; = Eo, eg = — Ej, e3 = Es is such an orthonormal frame. 


Hereafter, we suppose that e3 4 +E, = +2. In this case, we take e} = £5 


|e3 x Es] 
which is horizontal since (e,,e3) = 0, and obtain a natural orthonormal frame 
{e,, €2 =e3 X €1, e3} on M? x R which can be expressed as 


3 
(12) e:= > al E;,i = 1,2,3, (a) € SO(3). 
j=l 


By the choices e;, we have 
(13) a? = (e1, E3) = 0, a3 4 +1, (and hence) a3 Æ 0, 
and hence e; = atE, + af Ea with (at)? + (a?)? = 1. 
Furthermore, we can check that 
(14) fi=0, Vae = 0. 
In fact, a straightforward computation using (4) gives 
(15) 


3 . 3 . . 3 E . 
—fi 0 aE; = -fie = Verer = Va (2 aki) = 2 eilai) Eit Do aja, V g Ki. 
i=l i=1 


i=l i j=1 


Me 
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Using (2), the fact that a? = 0, and comparing the coefficient of E3 on both sides 
of (15) yields —f,a3 = e;(a?) = 0. From this we obtain fı = 0 since a3 4 0, and 
hence Vee = 0. 

It is not difficult to check that 
(16) E = a E: ae a; Es, E> = a; by + a? E, Es = Es 


defines another orthonormal frame on U x R C M? x R which satisfies 
[E1, Es] = f3E2, [E2, Es] = —fsEi, [E1, Ea] T, fE, 7 eae 

(17) Var = — fi, Via, E2 = fE, Vabi = — fE», Va, E2 = fF, 
all other Vg, E; = 0, i, j = 1,2,3, 


with 
(18) f = E,(aq) + Fx(aq) + fat, fs = af E3(aj) — aj E3(aj), 
and the only possible nonzero components of the (Ricci) curvatures given by 
Ric (£1, Ey) = Ric(B2, E2) = Rin = g(R(F1, E2) Eo, E1) = —E2(f) — P’, 
which is the Gauss curvature of M? on U denoted by KY? = —E,(f) — f?. 
We can also check that by introducing the new variables 0, a so that 
(19) a! = cos 6, a? = sin 0, a3 = sina, a3 = cosa, 


we can rewrite (16) and (12), respectively, as 


(20) E = sin 0E; — cos 0E, Ey = cos 0E; + sin 0E, E; = Es, 


e1 = cos 0E + sind Eo, 


(21) e = — cos q (sin OF, — cos 0 E2) + sin a F3, 
e3 = sin q(sin 0E; — cos 0 E2) + cos a E3, 
where 
at = cos, a? = sinf, a = 0, 
(22) ay = — cosa sin 0, az = cosa cos 0, a3 = sin a, 


az = sina sin ĝ, a} = — sina cos 0, a3 = cosa. 


It is also clear that the relationship between the orthonormal frames {e;} and 
{E;} is given by 


(23) e1 = Eo, eg = — cos aE; + sina E3, e3 = sin a E1 + cos a E3. 
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To compute the integrability data of the frame {e1, e2, e3}, we first use (2), 
(4), a? = fı = 0, and a computation similar to those used to obtain (15), to have 


€1(a3) = —oa3, €1(a3) = oa}, e2(a3) = 0, e2(a3) = 0, 


) 
(24) e3(a3) = —K2a3, €3(a3) = K203, K103 = (0 — f3)a3, fra, = 043, 


e1(aj) = —ajai f, e1(az) = (ai) f. 
Using (2), (23), (17), the 1st and the 2nd equations of (24) we have 


(25) fse2 + K1e3 = [e1, e3] = [e1, sin aE; + cos aE] 
= (o + f sin a cos a — cos? a E3 (0) Jez + (— f sin? a + sin a cos a E3 (0))eg. 
Comparing coefficients on both sides of (25) yields 
(26) fs =c + fsinacosa — cos? aE3(0), kı = —f sin? a + sin a cos a E3(0). 


Similarly, by computing 
(27) 
fze2 — 20e3 = [e1, e2] = [e1, — cos aF; + sin a E3] 
= (—f cos? a — sin a cos a E3 (0) Jez + (—o + f sin a cos a + sin? vE3(8))es, 


and comparing coefficients on both sides of this equation, we have 

(28) o = —f sina cosa — sin? aE3(0), fz = — f cos? a — sin a cos a E3(0). 
Using the 1st equation of (26) and the 1st equation of (28), we obtain 

(29) fs = —E3(0) = sin 0E3(cos@) — cos 0 E;(sin 0) = a? E; (at) — al Ez (a2). 
Using the 5th and 6th equation of (24), a simple computation yields 

(30) K2 = sin qez (cos a) — cos ae3(sina) = —e3(a). 


From (14), (26), (28), (29), and (30), we obtain the integrability data (10). 


Now we are ready to prove the following theorem which provides the main tool 
to prove our classification theorems. 


Theorem 2.6. Let n : M? x R > (N?,h) be a Riemannian submersion. Then, 
the orthonormal frame {e;} given by (9) is adapted to the Riemannian submersion 
m with the integrability data 

(31) 

fi =0, fo=—fcos*a, fg =0, o = —fsinacosa, kı = —f sin? a, k = —e3(a). 


Proof. By Proposition 2.5, we have V.,e; = 0. Thus, we can use Lemma 2.4 
to conclude that either V.,e2 Æ 0 in which case the chosen frame {e1, €2, e3} is 
adapted to the Riemannian submersion 7, or Ve e2 = 0. So, it suffices to prove 
that the orthonormal frame is also adapted in the the latter case: V.,e2 = 0, i.e., 
fo = 0. Using the assumptions (13), fo = 0, the 2nd and the 8th equations of (24) 
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we obtain ø = 0. Then, using the 4th equation of (5) with fı = f2 = o =a? = 0, 
we have either K™* = 0, or a3 = 0 and hence (a3)? = 1. 


For the case of a3 = 0 and (a3)? = 1, one applies the 7th equation of (24) 
with go = 0 to conclude that fs = 0, which implies that the chosen frame 
{e, = al E; + a} E2, e2, e3} is adapted to the Riemannian submersion. 


For the case K™” = 0 and a3 4 0,+1 (and hence a3 4 0,+1), we will prove 
that f3 = 0 and hence the orthonormal frame {e;} is adapted to 7. 

First, we note that the orthonormal frame {F;} given in (16) is a natural 
orthonormal frame with respect to the harmonic Riemannian submersion 


Tı : M? x R —> M°, ™(p,z) =p. 


By comparing the equations on the first line of (17) with those of (4) and using 
(29) and the second equation of (18), we find the generalized integrability data 
of {E;} to be 


(32) fi= f, =mi =R =F =0, fs = fs = —E;(0). 
Applying Lemma 2.4 with e) = E2,e. = Ei, we deduce that either fı = 


f #0 in which case the orthonormal frame {E;} is adapted to the Riemannian 
submersion 7, and hence fs = fs = 0, or fı = f = 0. For the latter case, we use 
the assumptions fı = f = 0, together with a3 4 0, +1 (and hence a3 Æ 0, 1), 
o = fz = 0 and (28) to have fs = —E3(0) = 0. Thus, in any case, the frame {e;} 
is adapted to the Riemannian submersion 7, and the generalized integrability 
data (10) reduces to the integrability data (31). 


Corollary 2.7. A biharmonic Riemannian submersion n : M? x R > (N?,h) 
with KM” = KN =0 has to be harmonic. 


Proof. By Theorem 2.6, the orthonormal frame {e;} is adapted to the Riemann- 
ian submersion 7 with the integrability data (31). Using the 4th equation of (5) 
and the assumption that KM = KN = 0, we have g = 0. From this and the 
4th equation of (31), we have either f =0 or sinacosa = 0. Therefore, we can 
obtain our corollary by the following cases: 


Case I: f = 0 and a3a3 = sinacosa # 0. Combining these, (31), and the 
assumption that K™ = KN =0 we have 
(33) 

at = fi = fo = fg =o = kı = K™ = K” =0, 03, a3 #0, +1, k3 = —ez(a). 
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Using (33) and (5) we obtain 
(34) €2(K2) = K2, e1(K2) = 0. 


Thus, in this case, if m is biharmonic with K = 0, then we have (4), (33), (34), 
and the biharmonic equation (8) reduces to 


(35) Ako = 0, 
and hence 
(36) €3€3(K2) = —K3. 


Applying ez to both sides of (36) we have 
(37) €2€3€3(K2) = —3K5. 


On the other hand, we use (4) and e2€3(K2) — €3€2(K2) = [e2, e3|(K2) = K2e3(K2) 
to obtain 


(38) €2€3(K2) = 3k 2€3(K2). 

By applying e3 to both sides of (38) and using (36) we get 
(39) €3€2€3(K2) = 3(€3(K2))” — 3K5. 

A further computation using (37) and (39) gives 


K2€3(€3(K2)) = [e2, e3](e3(K2)) 
= €2€3€3(K2) — €3€2€3(K2) = —3(€3(K2))?. 


This, together with (36), implies 
(41) 3(e3(K2))? = K3, 


Applying e3 to both sides of (41) and using (36) we get 


(40) 


(42) 10K3e3(K2) = 0, 


which implies either kọ = 0 or e3(K2) = 0. For the latter case, we use (41) to 
obtain k = 0. So in any case, the Riemannian submersion 7 is harmonic since 
Ki = kg = 0. 


Case II: a3 = cosa = 0 and a3 = sina = +1. In this case, (31) reduces to 


(43) fi = fo = fg =o =k. =0, kK = —f. 


Using these and (5) we have 


(44) €1(K1) = K?, €o(K1) = 0. 
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It follows that in this case, we can use (4), (43), (44) and K” = 0 to conclude 
that the biharmonic equation (8) reduces to 


(45) Akı =i): 


We make a computation similar to those used to compute (36)-(42) in Case I to 
have kı = 0. So a is harmonic since kı = ko = Q. 


Case III: ał = cosa = +1 and sina = 0. In this case, using the 5th and 6th 
equations of (31) we have «kı = «2 = 0 implying that the Riemannian submersion 
m is harmonic. 


Summarizing all the above cases we obtain the corollary. 


Now we are ready to give a first characterization of proper biharmonic Rie- 
mannian submersions from M? x R. 


Theorem 2.8. Let 7: M? x R > (N?,h) be a proper biharmonic Riemannian 
submersion. Then we have 

(a) The target surface is flat, and the adapted frame (9) has the integrability data 
fi = fo=ko =0 =0, kı =—f £0 solving the following PDE 


(46) Ak; = 0, i.e, AM’ f = 


where A“ denotes the Laplacian of M? and f = — sin @E,(0)+cos 0E2(0)+f sin 0 
is a function on U C M?, or, 

(b) The target surface is non-flat, and the adapted frame (9) has the integrability 
data 


(47) fi = ko =0, fo=—fcos*a, kı = —f sin?a, o = —f sinacosa 


satisfying fakio #0, K“ = e(f2) — f2 #0, o = —ei(a), elfa) = elki) = 
eilo) =e(K™’) = elf) = ela) = 0, i= 2,3, and kı solving the following PDE 


(48) Aki —Ki{-K + fz} = (0, 
which is equivalent to 
PL 


(49) a!” sin a cos? a + cos a(sin? a + 3)a’a”" + sin a(2 cos? a + 3)a” = 0, 


where a’, a” and a!” denote the first, the second and the third derivative of a 
along the vector field e1. 
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Proof. By Theorem 2.6, the local orthonormal frame given by (9) is adapted to 
the Riemannian submersion m with e3 being vertical and the integrability data 
given by 


(50) i = 0, fə = —fcos?a, fs = —E3(0) = 0, o = —f sin a cos q, 


kı = — f sin? a, Kp = sin ae3(cos.a) — cos qez(sin a) = —e3(a). 


We may assume that cosa = ał # +1, for otherwise e3 = +E3, m is the 
projection M? x R > M? followed by a Riemannian covering map and hence 7 is 
harmonic. The rest of the proof will be done by considering the following three 
cases. 

Case I: a3 = cosa = 0. In this case, it follows from (50) that fı = fo = k2 = 
o = 0, kı = —f #0. From these and (6), we conclude that N? has Gauss 
curvature KN = 0 and hence it is a flat surface. Therefore, biharmonic equation 
(8) turns into 


Ak, = 0. 


It is easy to see that the above equation is equivalent to A” : f = 0 since E30 = 0 
and hence Ez f = 0 and Af = AY’ f + E3E3f = AY f. 


Case II: a3 = cosa Æ 0, +1 and fz = 0. In this case, the second equation of (50) 
implies that f = 0, and hence we have fı = fo = fs = kı = 0 = KY? = KN 
0. Thus, we apply Corollary 2.7 to conclude that the biharmonic Riemannian 
submersion is harmonic in this case. 

Case III: a3 = cosa # 0,+1 and fə 4 0. In this case, the hypotheses can be 
summarized as: 


(51) a3, a} # 0,1, f2 #0, ay = fi = fs =0, e3(fi) = e3( f2) = 0. 


Claim 1: In this case, we have kọ = 0, eo(f) = e3(f) = ea(fo) = e3(fo) = 
€o(K1) = €3(41) = e(o) = e3() + e(K™) = e3( KM") = ez(a) = ez(a) = 
0, kac #0, and a’! = —o = fsinacosa, where a’ denotes the first derivative of 
a along e1 = Fy. 

Proof of Claim 1: Firstly, we use (51), the 7th and the 8th equation of (24) to 
have kia # 0. 

Secondly, by using (50) and the last equation of (51), we have e3(— cos? af) = 
e3( f2) = 0 and hence 


(52) e3(#1) = €3(— sin’ af) = e3[—(1 — cos” a) f] = —es( f). 
We use (50), the 5th and the 7th equation of (5) with a? = 0 to obtain 


(53) o? = ki fo, €2(K2) = K2, e2(0) = 2ko0. 
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On the other hand, note that a3 = sina, a3 = cosa, a straightforward computa- 
tion using the 6th equation of (24) gives 


0 = €3( f2) = e3(— cos? af) = —2Kq sin a cos af — cos? ae3(f). 
It follows that 
(54) e3(f) = —2kK2 sinacosaf __ 2K20 


cos? a ~~ cos? a’ 


Using (52), (54), the 1st equation of (53) and the 3rd equation of (5) with a? = 
fi = 0, we obtain 


e3(o2 e3(k ra 
e3(0) = fo) = heal kaf; 


(55) 


e1(k2) = e3(0) + K1K2 = Kof — k2 Sin? a f = k2 c08? a f = — ko f2. 
On the other hand, using ez(cosa) = e9(sina) = 0, f = -zZz and ez(0) = 
2k20 and a simple computation we have 
e2(f) = =e = aoe = ~~ kof, 
(56) €2( f2) = €2(— cos? af) = — COS ? aea( f) = = —2k2cos* af = 2kofo, 
€2(K1) = e2(—sin* af) = — sin? ae2(f) = —2K2sin? af = 2K1 ko. 


Note that K2e3(f2) = [e2,e3|(fo) = e2€3(fo) — e3e2( f2), a further computation 
using (56) gives 0 = —e3(2kef2) = —foe3(K2), which, together with fo 4 0, 
implies that 


(57) €3(K2) = 0. 


Applying the 2nd equation of biharmonic equations (8) and using (4), (53), (55), 
(56), (57) and (6) with fı = f = 0, we have 
(58) 
0 = -Ak + 2fre2(K1) + K1€2(f2) — Kika fe + Ko{—K% + f3} 
3 
= J {-exe;(K2) + Ve,e:(K2)} + 4k Ko fo + 2k Ko fo — Kika fo + Ko{—e1(fo) + 2f3} 
i=l 
= —K3 + 4k) Ko fo = Ko(—K2 + 4k fo) = Ko(—K2 + 407). 
Solving (58) we must have kz = 0 and «3 4 40°. Indeed, if K? = 40°, then using 
(57), we have 4e3(07) = e3(k3) = 0 and hence e3(c) = 0. But since f 4 0 and 
(55), we further have k2 = 0 and hence 40? = «3 = 0, a contradiction, since 
o #0. Therefore, a straightforward computation using (52)—(56), the 5th, the 
6th equation equation of (24) gives 


(59) 

eo(f) = e3(f) = = €3(cos@) = €2(K1) = e3(K1) = e3(0) = e2(0) = e2( f2) = 0. 
Since KY? = -Bf — P e1(f) — f2, we see that e3(K™”) = e.(K™”) = 0. 
From (24) and (59), we oe ole ) = e3(a) = 0, and a direct computation gives 
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osina = e;(cosa) = —sinaa’, i.e., a’ = —o = sinacosaf, where a’ denotes 
the first derivative of a along e; = Ey. For kə = fı = 0, biharmonic equation (8) 
turns into (48). A further computation using Claim 1 and (50), we see that (48) 
turns into (49). 


Finally, we show that (N?, h) is not flat. If otherwise, K” = 0, i.e., e1(f2)—f? = 


0, using the 4th equation of (5) with o = — f sina cosa, we immediately have 
(60) KM = -3P sin? a. 

On the other hand, we know that M? has Gauss curvature 

(61) KM = —By(f) — P =—e(f) — P. 

Using Equations (60) and (61), we obtain 

(62) f =3f? sin’ a — f’, 


where f’ denotes the first derivative of f along ey. 
Substituting a’ = sinacosaf (i.e., f = ) into (62) and simplifying the 
resulting equation we get 


(63) a" cosa — a” sina = 0. 
One applies e; to both sides of (63) and simplifies the resulting equation to have 


3 


(64) a” cosa — 3a'a” sina — a” cosa = 0. 


Adding (49) to a (— sin a cos a) multiple of (64) and simplifying the results with 
a’ Æ 0 yields 


(65) a" cosa(4sin? a + 3) + a” sin a(3 cos? a + 3) = 0. 


Similarly, adding a (—4sin? æa — 3) multiple of (63) to (66) and simplifying the 
results gives 


(66) a° sina(9 + sin? a) = 0, 


it follows that a is a constant, a contradiction. Then, we must have KN Æ 0. 


Summarizing all results in the above cases we obtain the theorem. 


We now give a characterization of proper biharmonic Riemannian submersions 


from M? x R by using the local coordinates as follows 


Theorem 2.9. Jf 7: M? x R > (N?,h) is a proper biharmonic Riemannian 
submersion from the product space, then 
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(i) The target surface is flat, and locally, up to an isometry of the domain and/or 
codomain, m is the projection of the special twisted product 


(67) T: (R, e2P(@9) dr? + dy? + dz”) > (R’, dy? + dz”), x(x,y,z) = (y,z), 


with py # 0 being a harmonic function on (M?, e??@)dx? + dy”), i.e., it solves 
the PDE 


(68) ADy := Pyyy + PyyPy + EPE (Dewy — PryPr) = 0. 


or, 
(ii) The target surface is non-flat, and locally, up to an isometry of the domain 
and/or codomain, the map can be expressed as 


m: (RÌ, 2?) dx? + dy? + dz?) > (R°, dy? + eA d¢?), 
T(£,Y, z) = (y, F (z z fer@dz)), 


where p(x,y) = In|tana(y)| + y(x), A = In|sina(y)| + w(¢) with the functions 
plz), w(d) and nonconstant function F(u) satisfying F(z — f e? dr) = e~ 
and z — f e@dr = f edd, and a(y) is the angle between the fibers of x and 
E; = 0, solving the ODE (49). 


(69) 


Proof. First of all, note that by Theorem 2.8 the local orthonormal frame {e; = 
Ey, e2 = — cos aF +sina E, e3 = sinaE,+cos aE} is an adapted frame of the 
Riemannian submersion 7, and that the vector field e} = E> is a geodesic vector 
field on M?. It is well known that we can choose local semi geodesic coordinates 
(x,y) on M? so that Ey = e~P@9,, Ey = Oy, and the metric on M? takes the 
form 7?) dx? + dy?. It follows that the product manifold M? x R can be locally 
represented as (U x R C M? x R, €??@) da? + dy? + dz?). 

For Statement (i), since the target surface is flat, it corresponds to the Case 
I in the proof of Theorem 2.8, i.e., the frame {e, = E, e2 = Ez, e3 = E} is 
an adapted frame to the Riemannian submersion m with the integrability data 
{fi = fo= kg =o =0, kı = —f #0}. Note that, by o = 0 and the 4th equation 
of (5), the horizontal distribution of the Riemannian submersion is integrable with 
flat integral submanifolds. So locally, up to an isometry of the domain and/or 


the target manifold, the Riemannian submersion is the projection along the fibers 
(i.e., the integral curves of E, = e~?!0, to the integral submanifold, and hence 
can be described by (67)). It is easily checked that in this case kı = —f = —p, 
and Equation (46) reduces to (68). 

For Statement (ii), Theorem 2.8 implies that in this case, the target surface 
is non-flat, cosa Æ constant depending only on variable y, and {e; = Es, e2 = 
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—cosaFk, + sinak3, e3 = sina, + cos aF} is an adapted frame of the Rie- 
mannian submersion 7 with integrability data (47). 

Note that the vector field e4} = E> = Oy is a basic vector field to the Riemannian 
submersion 7. It is well known that there is a local vector field £1 on (N?, h) whose 
integral curves are geodesics on (N, h) such that dr(e1) = £1. It follows that we 
can choose a local semi geodesic coordinates (y,¢) on N so that metric takes 
the form h = dy? + e?%-*)d¢?, and the orthonormal frame c; = 0,,€2 = e~*04 
satisfying dm(e,) = €1,d(e) = £2 = € 04. 

Summarizing the above, we conclude that in this case, up to and isometry of 
the domain and/or target manifold, the Riemannian submersion can be expressed 
as 
(70) 

m: (RË, cP) da? + dy? + dz”) > (R°, h = dy? +e?) dg’), (x,y, z) = (y, Q), 
where @ = (x,y, z) is a function to be determined. 


Now we are to determine the functions ¢ = (x,y,z), p(x, y), and A(y, Q). 
A straightforward computation gives 


(71) Fy, ome, + Fo oTe = [eE1, E2] = —AyE2, 


which implies 
Fy Fa 0, Fy = — Ày. 


Hence, we have 


(72) Ay = —Fhon = — f2 = f co? a = raou 
Integrating both sides with respect to y yields 
(73) Ay, d) = S SRA dy + w(4) = In|sina(y)| + w(¢), 


where w = w(@) is an arbitrary function on @. 


Noting that p, = f(y) = 2(Y) _ does not depend on x we have 


sin @ cos Q 


(74) p(a,y) = J fy)dy + p(z) = f <2? dy + g(x) = In| tan a(y)| + p(z), 


where y = ¢(z) is a function on z. 


To determine the component function (x,y, z), we use e1 = Oy, e2 = — cos ae PPO, + 
sin adz,e3 = sinae?)0, + cosad,, and dr(e1) = €1 = Oy, dr(e2) = €2 = 
e 04, and d7(e3) = 0 to have 
Oy = €4 = dr(e1) = Oy + Py Op, 
(75) ae =£, = dr(ez) = (—cosae Pd, + sin ah, )Og, 
0= dr(e3) = (sinae?*™) p, + cos ag, )Oz 
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By comparing both sides of the 1st equation of (75), one finds that ¢, = 7 =), 
which means that the function ¢@ does not depend on y, i.e., ¢ = ¢(2, z). 
Comparing coefficients of both sides of the 2nd and the 3rd equation of (75) 
separately, we get 


(76) —cosae Pg, + sinad, =e, sinae? ¢, + cosad, = 0. 


Recall that p(z,y) = In|tana(y)| + y(x), A = In|sina(y)| + w(¢), and the 
fact that ¢(x,z) is nonconstant since da(e2) 4 0, we use the method of the first 
integral to solve the 2nd PDE of (76) to have 


(77) lx, z) =F (: = / tr) : 


where F = F(u) is a nonconstant differentiable function. Substituting this into 
the lst PDE of (76) we have F'(z — fe?™dzx) = e~” . It follows from this 
and (77) that dọ = dF = F'du = e~’ du and hence e” do = du implying that 
u =z — fe?@dxr = fedo. This completes the proof of Statement (ii). 


Applying Theorem 2.9, we immediately have the following corollary which char- 
acterizes a proper biharmonic Riemannian submersions from a product manifold 
onto a non-flat surface as a special map determined up to an arbitrary function 
between two special warped product manifolds with the warping functions solving 
an ODE. 


Corollary 2.10. A proper biharmonic Riemannian submersion m : M? x R > 
(N?,h) from product manifold into a non-flat surface is locally, up to an isometry 
of the domain and/or codomain, n is a map between two special warped product 
spaces given by 

: (R®, tan? a(y) dt? + dy? + dz?) > (R?, dy? + sin? a(y) dy?) 

m(t,y,z) = (y,z — t), 


T: 


~ 


(78) 


where a(y) is the angle between the fibers of n and E3 = 2 solving the ODE 


(79) a” sin a cos? a + cos a(sin? a + 3)a'a” + sin a(2 cos? a + 3)a” = 0, 
Proof. This follows from Statement (ii) of Theorem 2.9 and the coordinate changes 
t = fedr, y = y, z = z in the domain and y = y, Y = fe’ dd in the 
codomain. 


Remark 3. (A) Note that it follows form [1] (Corollary 3.2) that the Riemannian 
submersion given by the projection of the twisted product m : (RÌ, e??(%7)da? + 
dy? + dz”) > (R?,dy? + dz?), a(x,y,z) = (y,z) is biharmonic if and only if 
Ak, = 0, Akg = 0. In the case of (i) in Theorem 2.9, these reduce exactly to (68). 
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Thus, Statement (i) of Theorem 2.9 not just characterizes proper biharmonic 
Riemannian submersions from product manifold to a flat surface locally but also 
recovers a special result in [1] (Corollary 3.2). 

(B) The biharmonicity of the Riemannian submersion defined by the projection 
of warped product r : (RÌ, e??7) da? + dy? + dz?) > (R°, dy? + dz”), m(x, y, z) = 
(y, z) had been studied in [2, 15, 21, 12, 1]. 


(D) We would like to point out that, unlike in case (ii) in Theorem 2.9 and 
Corollary 2.10 where the function p(z, y) can be proved to be independent of x 
variable so the metric is of a warped product type, in case (i) the map is complete- 
ly determined (to be the projection) but we do not know whether the function 
p(x, y) is independent of x variable or not. 


When the target surface is flat, it is easy to have many examples proper bi- 


harmonic Riemannian submersions 7 : M? x R — R? from the projection of 


the warped product spaces see e.g., [2, 15, 12, 1]. For example, the following 
projections are proper biharmonic Riemannian submersions: 


(i) 7 : (R? x R, (cosh y)'de? + dy? + d2?) > (RÈ, dy? + dz”), (x,y, 2) =(y,2), 
(ii) 7: (R? x R, y*dx? + dy? + dz”) > (RÊ, dy? + dz), n(x, y, z) = (y, z). 


When the target surface is non-flat, we would like to point out that there exist 


many local proper biharmonic Riemannian submersions 7 : M? x R > (N?,h). 


H 
In fact, by introducing new variable u(a) = e. we have a'a” = ua’, al” = 


(u’ + 2u?)a’? and hence (79) reduces to a Riccati equation as 


1 | 2 | sin?a+3 | 2cos*at3 _ 
(80) u (a) F Que 4 sinacosa’ T ~ cos? a = 0, 


any solution of which gives a family of locally defined proper biharmonic Rie- 
mannian submersions M? x R > (N?,h). 

Finally, note that it was proved in [21, 22] that a proper biharmonic Riemannian 
submersion M?(c) xR exists only in the case when c < 0, and m : (RÌ, e2¥~“¥da? + 
dy? +dz?) > (R°, dy?+ dz”), m(x,y, z) = (y, z) is an example. Now we can prove 
that, up to isometry, this is the only one. 


Proposition 2.11. A Riemannian submersion n : M?(c) xR — (N?,h) is proper 
biharmonic if and only ifc < 0, (N?,h) is flat, and, up to an isometry, the map 
can be expressed as m : H?(c) x R > R? with m : (R°, e2Y~°4dx? + dy? + dz?) > 
R?, dy? + dz”), n(x,y,z) = (y, z). 


“—~ 


Proof. Firstly, it follows from [22] that a proper biharmonic Riemannian sub- 


mersion m : M?(c) x R > (N?, h) from a product space exists only in the case: 
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H?(c) x R > R? with c < 0. 


Secondly, by Theorem 2.8 and 2.9, we know that locally, up to an isometry 
of the domain and/or codomain, a proper biharmonic Riemannian submersion 
m : H?(c) x R > R? with c < 0 is expressed as 


(81) T: (R, erv) dr? + dy? + dz”) > (R’, dy? + dz”), x(x,y,z) = (y,z), 


and the orthonormal frame {e; = Oy, e2 = Oz, e3 = e ?0,} is adapted to the 
Riemannian submersion 7 with the integrability data fı = fo = k2 = o =0, kı = 
—py #0. It is easily checked that in this case, (5) reduces to 


(82) ei(kı) = k? +c, €2(K1) = 0, 
and biharmonic equation (8) reads 

(83) Akı =0. 

A straightforward computation gives 


AK, = e1e1(K1) + €3€3(K1) = Veels) = Vex€2(K1) = Veg€3(K1) 
= e (K? +c) + €3€3(K1) — K1e1(K1) = €3€3(K 1) + K3 + CK. 


Substituting (84) into (83), we have 


(84) 


(85) €3€3(K1) = —K3 — CK}. 


Applying e3 to both sides of the 1st equation of (82) and using the fact that 
€1€3(K1) = [e1, e3|(K1) + e3€1(K1), we have 


(86) €1€3(K1) = 3K1€3(K1). 
Using (82), (85), (86), and a direct computation we get 
(87) 
e1e3{e3(K1)} — e3e1{e3(K1)} = [e1, e3 {e3(K1)} = K1e3¢3(K1) = —K} — cr?, 
and 
(88) e1€3{€3(K1)} — e3e1{e3(K1) } = er {ess (ei) } — es{eres(K1)} 


= —cki — 4c — 3€3(K1). 
Comparing (87) with (88), we get 
(89) 3e3(K1) = KÍ — 4c. 


Applying ez to both sides of (89) and using (85) to simplify the resulting equation 
we have 


(90) K1(5K? + 3c)e3(K1) = 0, 
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which implies e3(K;) = 0. Substituting this into (85) and using that fact that 


CRY 


# 0 we obtain 


(91) k? = —c > 0, 


which implies that 


(92) p? = —c, (and hence) p, = +y—c. 
It follows that 


(93) 


plz, y) = £V—cy + p(z), 


where (x) is a an arbitrary function. 


So we conclude that up to an isometry of the domain and/or codomain, a 
proper biharmonic Riemannian submersion m : H?(c) x R > R° with c < 0 is 
expressed as 


(94) 


Thus, we obtain the proposition. 


m: (RS, e?e tdg? + dy? + dz”) > (R°, dy? + dz”), n(x,y, 2) = (y, 2). 


Remark 4. By Proposition 2.11, for c > 0, there exists no proper biharmonic 
Riemannian submersion 7 : M?(c) x R > (N?,h) no matter what (N?,h) is. 
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